Protected quantum computation with multiple resonators 
in ultrastrong coupling circuit QED 
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We investigate theoretically the dynamical behavior of a qubit obtained with the two ground 
eigenstates of an ultrastrong coupling circuit-QED system consisting of a finite number of Josephson 
fluxonium atoms inductively coupled to a transmission line resonator. We show an universal set 
of quantum gates by using multiple transmission line resonators (each resonator represents a single 
qubit). We discuss the intrinsic 'anisotropic' nature of noise sources for fluxonium artificial atoms. 
Through a master equation treatment with colored noise and manylevel dynamics, we prove that, 
for a general class of anisotropic noise sources, the coherence time of the qubit and the fidelity of 
the quantum operations can be dramatically improved in an optimal regime of ultrastrong coupling, 
where the ground state is an entangled photonic 'cat' state. 
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The study of quantum dccoherence is believed to be 
crucial in order to understand the transition from the 
microscopic quantum world to the macroscopic classical 
one. Moreover, a control and limitation of decoherence 
is essential towards the realization of a robust, scalable 
quantum computer. The study of cavity QED systems 
in atomic physics [l[ has led to spectacular fundamen- 
tal investigations of non-unitary evolution due to deco- 
herence mechanisms. In particular, it has been possible 
to observe the fragility of states of the form \^ ca t) = 
■^{\a) P hot\g)at + \ae l11 ) pho t\e)at} (usually dubbed 'cat 
states'dl), where \a) v hot is a coherent photon state with 
a large mean photon number \a\ 2 ^> 1, \ae lr, ) v hot is an- 
other coherent state with a phase difference rj, while \g) a t 
QeVnt) is the ground (excited) state of a two- level atom 
|2|, |3| . These states have been prepared in a cavity QED 
system well described by the Jaynes-Cummings model, 
where the ground state is \0} p hot\g}at, i-e., the vacuum of 
photons times the atomic ground state. Recently, a grow- 
ing interest has been generated by the so-called ultra- 
strong coupling regime of cavity (circuit) QED both the- 
oretically [1, [SB-El and experimentally (li - flit . Such 
a regime is achieved when the vacuum Rabi frequency 
r^o, which quantifies the coupling between one photon 
and one elementary matter excitation, is comparable or 
larger than the cavity (resonator) photon frequency uj cav 
. In such a regime, the Jaynes-Cummings model based 
on the rotating wave-approximation (valid for small ratio 
flQ/Lo cav ) breaks down. In particular, the ground state of 
the system is no longer the standard vacuum: recently, it 
was shown [Toj j that in the limit of very large coupling 
the ground state can become quasi-degenerate with the 
entangled structure: 

|tt G ) ~ (\a) ph nf =1 \ +)j + (_i)"| - a ) ph uf =1 \-) 3 ) 

(1) 

where N is the number of atoms embedded in the cavity 



resonator, \a) p h is a coherent state for the photonic field 
which satisfies \a\ ~ ^/NQo/uj cav , and |±)j are pseudo- 
spin polarized states for the j-th artificial atom, which 
are defined in the following. For each two-level system 
{|e)j,|<7)j} one can introduce the Pauli operators a 3 x = 



\e)j(9\j + \g)j( e 
(2|e> 3 -(e 



*(Is>3 



and &{ 



1). With a light-matter coupling Hamiltonian 



of the form H, 



coupling 



-EL A, 



h.c, (Xj being 



the local coupling strength and a the photonic bosonic 



annihilation operator), |±)j 
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\e)j ± \g)j) are the 



eigenstates of a 3 x . Interestingly, the (orthogonal) first 
excited state has the similar form: 

\* E ) ~ i (|a) ph njL 1 |+) j - (-in - a) ph nf =1 \-),) 

(2) 

In this letter, we show how \^g) and \^e) surprisingly 
can form a robust qubit, whose decoherence can diminish 
while increasing the 'size' of the corresponding photonic 
'cat' states (see Fig. 1). Moreover, we also provide a 
universal set of quantum computation gates and demon- 
strate via a thorough master equation treatment the fi- 
delity enhancement in a regime of ultrastrong coupling. 



The energy difference <5 between the two considered states 
diminishes exponentially 0, ITo| with the vacuum Rabi 

coupling, namely S ~ u eg exp(— 2—^-N), where u> eg is 
the frequency of the single atom two-level transition (set 
to be equal to the cavity mode frequency). Either in 
the ultrastrong (fl /ui e g — > +oo) or 'thermodynamic' 
(N — > +oo) limit, the two states become degenerate. In 
the ultrastrong coupling limit, the other excited states 
are much higher in energy, separated by a frequency gap 
A ~ Lu eg >> S. Importantly, these interesting features 
can not be obtained in every ultrastrong light-matter 
coupled system. In particular, using the Pauli matrix 



2 



/I 




1 Ultrastrong circuit-QED Resonator =1 'vacuum' qubit 

= = ^{|a> | + +) - 1 - q) | )} 

^1*> = ^{I«) | + +> + l -«>!-->} 




FIG. 1: Description of the considered system. The building 
block is a superconducting transmission line resonator em- 
bedding N Josephson atoms (N — 2 in the sketch here). 
By choosing judiciously the type of artificial atom (the de- 
picted circuit represents fiuxonium atoms inductively coupled 
to the resonator) , the first two ground levels of the resonator 
are entangled states (\a) is a photon coherent state, |±) is a 
Josephson junction state 'polarized' along the pseudospin x- 
direction). One resonator represents a single qubit: a register 
of M qubits is given by M resonators. 



language, the 'direction' of the bare atomic Hamiltonian 
must be orthogonal to the one of the light-matter inter- 
action Hamiltonian [l^. This is the case in the following 
spin-boson Hamiltonian: 



H/h = uj cav a)a + 
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N 



N o 
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(3) 



In the present letter, we limit our description to a single 
bososic mode and an uniform light-matter coupling , but 
all the following results may be generalized to several and 
spatially non uniform modes 

It has also been shown recently that in the ultrastrong 
coupling regime, the quasi-degeneracy of the states |\?g) 
and \*&e) is robust with respect to a local and static per- 
turbation of the type = J2j=i hy-j&y.j + hz,j&z,j 
where h y _j and h z j arc random perturbation amplitudes 
Q. The reason is that in the subspace |\&e)} 
such perturbation couples (at the iV'^order) coherent 
states of opposite phase | — a) and \a). The effect of 
the perturbation is proportional to the overlap (— a\a) = 
cxp(— 2|a| 2 ) ~ cxp(— 2—^-N). Indeed, the stronger is 
the coupling D,q or the larger is the number of artificial 
atoms TV, the larger is \a\ 2 , the 'size' of the photonic 'cat' 
states \^fg) and \^e)- Importantly, the protection is not 
complete [19(, because these states are not robust with 
respect to noise terms like fff ert = Ylj=i h>x,j&x,j and 
H? ert = h a a + h*a\ namely the noise in the direction of 
the light-matter coupling and the noise associated to the 
resonator field. However, if in a superconducting system, 
perturbations like -fff e 2 ' ' happen to be the dominant ones, 



the lifetime and the fidelity of the quantum operation in- 
volving the states \^g) and \^e) can be dramatically 
improved by increasing il /uj eg and/or N. 

In fact, among the different flux Josephson atoms 17, 
18, 2(J [2l||, this noise anisotropy appears to be realistic 
at least for a fiuxonium (20j . Under conditions detailled 
in |2fJ, its Hamiltonian can be written as: 



Hp 



AE Cj Nj 



Er 



Ejcos{(pj + § ext ) (4) 



The Hamiltonian parameters are subject to noise fluctu- 
ations : $>ext = 7T + A§ ext with A$> ext some flux noise 
(in units of $ = h/2e), Ej = Ej + AEj with AEj = 
A/o/$o proportional to the critical current fluctuation, 
Nj = Nj + AN Q , ANq being the charge offset fluctua- 
tion. One can also introduce some capacitive and induc- 
tive noise E Cj = E Cj + AE Cj and E L j = E L j + AE Lj . 
When the fluctuation sources are off, the first two eigen- 
states of the fiuxonium are very well isolated from the 
higher states provided that Ej ^> Elj and Ej ^> Ecj- 
Then, the Hamiltonian (|4|) reads Hp — h(uj e g/2)a z in 
the basis of the two first eigenstates which are symmet- 
ric and antisymmetric superpositions of clockwise and 
anticlockwise persistent current states. On the same ba- 
sis tpj ~ -(foi&x and Nj ~ £ff^ipoio- y (where ip i -i 3). 
The fluctuations produce (at the first order) the pertur- 
bation: 

H F ,pert/h— A§ ex tSin((poi){Ej/h)a x + ANoifiOlUegVy (5) 

, s dw eg AI 



duj eg AE Cj 
-hjc 



dE Cj 



E, 



Cj 



duj eg AE L . 
y -E L ^^)a z 



dE L 



Ei 



The spectral density of the flux noise is typically 

S A^ t ~ W~ 6 /VHz [H,!!!- The critical current noise 
AIq/Iq = AEj/Ej, which is also believed to follow a 



1// law [24|, |25|, has been recently measured [26| in a 



1/2 



fiuxonium : S 



J AEj/E, ~ 3.10~ 5 /\/-ffz. It proves that the 
dissipation due to the a z channel is much larger than 
the a x channel contribution. To study the behavior of 
the qubit {I^q), in the presence of dissipation, we 

used the master equation 27J : 



^ = 1 [H , p] + Urf>Sr + S r pUl - S r U r p - pU}S r 

r=r v ,r f 

N 

+ ^ ^ UmpSm + S m pU} n - S rn U rn p - pU^Sm (6) 

j=l m=x j ,y j ,z j 

where p is the density matrix, H refers to Hamiltonian ((3]) 
and where the 'jump' operators are S rv = a + cv, S r , = 
i(a— at), S Xj = S Vj = a y , S Zj = a 3 z . Moreover [28), 



Uk= / Mr)e-^ HT S k ei HT dT,(7) 
Jo 

r k (u){n k (u;)e iUT + [n k (u) + l]e- iUT }aw, 



3 



for k = r v ,rf or k 



Here we consider the zero temperature limit [2 9f. where 
the spectral functions Tk(w) must vanish for u < be- 
cause they are proportional to the density of states (of 
the baths) at energy hu. For sake of simplicity, we have 
set Tfc(w) = Tfe for u € [0;w c ] and Tk(w) = elsewhere 
Vfc, with w c an upper cut-off which is consistent with de- 
creasing spectral noise. Finally, one must include many 
excited states in the master equation treatment. To in- 
vestigate the robustness of the coherence between the 2 
quasi-degenerate vacua \&g) an d |*e) , we have studied 
the non-unitary dynamics of the initially prepared pure 
state |* ) = cos(0)|*.e) + sin(6')e^|* G ) in presence of 
anisotropic Josephson dissipation rates T y ,T z 3> T x and 

for several cavity loss rates T r /uj eg = T r Ju eg = T rf /u eg FIG. 2: Coherence time in units of l/uj eg calculated via the 




(see caption of Fig. 2). Our simulations plotted in Fig. 
[2] prove that the coherence time increases while increas- 
ing the normalized vacuum Rabi frequency £lo/uj eg . In- 
deed, if the dominant dissipation channels are along the 
y and z directions, their effect decreases as exp (— 2|a| 2 ) 
where a = y/NClo/^cav Hence, the coherence time 
is enhanced exponentially before reaching a saturation 
value given by r r , T x and eventually decreasing with the 
usual power law of cat states. The location of the coher- 
ence time peaks with respect to the photonic amplitude 
a = V^V^o/^eg is almost independent of the number of 
atoms N (see top right panel of Fig. 2), indicating that 
a is the relevant dimensionless parameter for the protec- 
tion. Depending on r r and (see bottom right panel 
of Fig. 2), the maximum coherent times have a different 
behavior versus N. For smaller values of F x , the protec- 
tion increases monotonically with N > 2 (we have been 
able to calculate up to TV = 5). For larger values of T x in- 
stead the maximum of the coherence time is achieved for 
N = 1. Finally, since the number of photons (n) of \^g) 

and \ty e) increases like r^P-N (see Fig. |2J, we conclude 
that there is a regime where the larger is the number 
of photons in \^g) and \^e), the stronger is their ro- 
bustness against decoherence contrary to the usual cavity 
QED 'cat' states obtained when Vt a /io eg << 1. In- 
deed, it is well known[H-l3l.[30| that the coherence time of 
those standard 'cat' states decreases monotonically while 
increasing their size. 



Now, we show how to obtain an universal set of 



master equation ((6]l for 



ujcav and with the initial state 



gates for quantum computation [3l| using the two states 
I^g) and \*&e) as computational basis for the qubit 
and we will study the fidelity of such quantum opera- 
tions. One begins by showing how to get the dynam- 
ical gate e~ 4eiSl in the basis \^g) and \^e), where 
= + |*e)(*g| is the x-direction Pauli ma- 

trix associated to this (collective) vacuum qubit. To do 
so, one can add a coupling between the flux of one Joseph- 
son atom embedded into the resonator (for instance the 



|*o) = cos(0)|*s) + sin(6')e i<# '|*G). Results are averaged 
over the possible initial values for and <j>. Left panel: co- 
herence time versus the normalized vacuum Rabi frequency 
Oo/weg for one atom (N = 1) with Josephson loss rates 
{F x ,r y ,r z } = tj e9 {10~ 6 ,10~ 3 ,lCr 3 }. The different cavity 
loss rates: r r /tv eg = lCT 6 ,lCr 7 ,0 correspond [lj to differ- 
ent quality factors Q = uj eg /(47rr r ) ~ 10 5 ,10 6 ,oo. Inset : 
the number of photons (n) = \a\ 2 — (a/a) is plotted versus 
Qo/uje S for N = 1. Top right panel: coherence time for N = 1, 
2 and 3 atoms for r r /cj e9 = 10~ 6 and with a lower anisotropy 
in the atomic loss rates: {T x ,F y ,r z } = w e9 {lCn 5 , 1CT 3 , 1(T 3 } 
as a function of the photonic amplitude a = y/NClo /oJ eg . Bot- 
tom right panel: maximum coherence time as a function of 
the number of atoms N for different values of Y x , hence for 
different noise anisotropy. 



first atom) and an external, classical and tunable mag- 
netic field & s (t). This leads to an additional Hamiltonian 
term of the type M$> s (t)0j = C(t)a x where is the flux 
across the Josephson junction of the first artificial atom. 
Such perturbation lifts the degeneracy of the fundamen- 
tal subspace so that the new two first eigenstates are 
|+)| + a) and |-)| -a) with a splitting S(t) = 2C(t) and 
where we have replaced H^ =1 \zL)j by |±) to simplify the 
notation. By adiabatically shaping the time-dependence 
of C(t) it is possible to create a dynamical gate e' 16 ^* 
with 9 X = J Q T C(t)dt with [0;T] the gate time interval. 

Now, we show how to get a second single-qubit gate, 
namely e - lS ^- where S z = 2\^ e )(^e\ - 1. |* G ) and 
\^e) have an energy splitting S exponentially decreasing 
as a function of f^o- By modulating in time Qq, one 
gets the desired quantum gate. Acting adiabatically, the 
rotation angle will be 9 Z = L T S(t)dt = S(flo(t))dt. 
Even without optimizing the temporal shape of t — > f^o (t) 
, excellent fidelities can be reached. For instance, for the 
Z-Pauli gate (corresponding to 8 Z = 7r/2), with one atom 
and for a linear back and forth between flo/^eg = 2 and 
fl /uj eg = 1.3, fidelity > 99.9% is obtained for a typical 
time T ~ 300/u; e£ , in presence of realistic dissipation. 
In practice, to modulate in situ r2o(t), one can use an 
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intermediate loop between the resonator and the artificial 
atom with a tunable magnetic flux through it 0, [l2| . 




FIG. 3: Left panel: fidelity of the single qubit X rotation 
gate (for 6 X = 7r/2) versus Q.o/u) eg . Master equation ((6]l was 
used with a time-dependent Hamiltonian, H(t) — H(t = 0) + 
C(t)a x with H(t = 0) the initial spin-boson Hamiltonian Q 
with N = 1 (black solid), N = 2 (blue dashed) and N = 3 
(red solid) Josephson atoms in the resonator. Inset: time 
evolution of C(t). Right panel: fidelities of the 2-qubit gate 
e -i6a 13 t ai igit X3 for _ 7r /2 with respect to the vacuum 
Rabi Frequency in each resonator in which there is 1 atom 
embedded. The 2-qubits coupling constant C 12 (t) follows the 
same time evolution as C(t) in the inset. Note that we have 
included noise in the mutual coupling, via the jump operator 

S X 12 



i&x 2 an d the loss rate T x i2 



In order to get a complete set of quantum oper- 
ations, one needs to perform a 2-qubit control gate. 
Here, we will describe how to obtain the conditional 



{{|* G )l, \* E )l} ® {|¥ G > 2) l**M} = {^(| 



a)i ± 



|-)| - a)i) ® ^ (|+)| + a) 2 ± |-)| - a) 2 )} where 1 (2) 
stands for the resonator number. For our goal, one way 
is provided by a direct magnetic mutual coupling [32j 
M 12 (£)<y3j<^, between 2 fluxonium atoms (one in each 
resonator), giving the Hamiltonian H12 = Hi + H2 + 
C 12 (t)al x a\ 2 , where Hi (H 2 ) stands for the spin-boson 
Hamiltonian §5§ for the resonator 1 (2), while a\. x (resp. 
cr\ 2 ) stands for the x-Pauli matrix acting on the first 
two levels system of the resonator 1 (2). Applying such 
a perturbation will partially lift the 4 times degener- 
acy of the fundamental subspace so that the two states 
(|+)|+a)i®|+>|+a>2 and |-)|-a)i<g>|-)|-tv) 2 ) will have 
a different energy than the states (|+)| + a)i ® | — )| — a) 2 
and |— ) I — a) 1 <E) |+) | +0)2)- Fidelity of that operation for 
8 Xl2 = tt/2 is given in the right panel of Fig. [3] in presence 



of dissipation, showing again the enhancement for in- 
creasing values of the normalized vacuum Rabi frequency. 
Other proposals for the practical coupling between the 2 



resonators could be envisaged [33l 13-41 ] . Concerning the 
read-out of our qubit, this can be done by a projective 
measurement on the states |+)| + a) and |— )| — a): the 
flux across the Josephson junctions is polarized and can 
be in principle measured via the surrounding quasi-static 
magnetic field. 

In conclusion, we have shown that it possible to 
considerably enhance the coherence times of a qubit 
given by the first two eigenstates of a circuit QED 
system in the ultrastrong coupling regime: such states 
are entangled states of photons and polarized Josephson 
atomic states, which are robust with respect to a general 
class of 'anisotropic' environment. In our proposal, the 
resonator is used to protect quantum information 35 , 3(| 
, contrary to the approach 37 -39| where it acts as a 
bus joining several embedded Josephson qubits. The 
present work shows that the qualitative modification 
of the quantum ground state in ultrastrong coupling 
circuit QED can have a significant impact on the deco- 
herence and manipulation of quantum states in multiple 
resonators. We would like to thank M.H.Devoret for a 
critical reading of the manuscript and useful discussions. 
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